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Model  number  LS-ll-DR-2  and  a set  of  triaxial  Schonstadt  600  milligauss 
magnetometers.  An  attitude  determination  problem  was  created  when  the 
telemetered  lunar  sensor  angular  signal  was  saturated  for  nearly  the  entire 
flight;  however,  the  times  of  the  lunar  readout  were  accurately  measured. 
Vehicle  procession  and  spin  l'.tes  were  determined  from  corrected  magnetometer 
data  together  with  the  available  lunar  information.  Using  the  above  attitude 
information  in  conjunction  with  corrected  magnetometer  data,  standard 
attitude  determination  techniques  discussed  in  [2,3]  were  modified  to 
generate  lunar  pitch  angles. 
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The  knowledge  of  the  ionosphere  has  been  greatly  increased  in  recent 
years  through  the  study  of  data  provided  from  rocket -borne  sensors.  For 
certain  sensors,  meaningful  interpretation  of  the  instrument’s  output  can 
be  accomplished  only  when  the  orientation  of  the  sensor  or  its  transporting 
vehicle  is  known.  The  orientation  of  the  sensing  axis  of  an  instrument  with 
respect  to  a fixed  coordinate  system  in  space  is  referred  to  as  its  attitude, 
and  from  the  attitude  of  a sensor  one  may  determine  the  angle  it  makes 
with  any  other  vector  such  as  the  rocket's  velocity  vector  or  the  Garth's 
magnetic  field  vector. 

The  first  section  of  this  report  applies  the  data  reduction  techniques 
described  in  [1]  to  obtain  the  attitude  and  requested  angular  outputs  for 
most  rockets  investigated.  These  rockets  provided  gyroscopic  outputs  which 
measure  the  vehicle's  roll,  pitch,  and  yaw  with  respect  to  a coordinate 
system  fixed  at  launch. 

The  analysis  and  data  reduction  techniques  discussed  in  the  remaining 
sections  of  this  report  were  utilized  to  determine  vehicle  attitude  and 
angles  of  attack  of  the  on-board  probes  of  the  Ute-Tomahawk  Rocket, 

A09. 209-1,  fired  16  April  1973  at  White  Sands  Missile  Range,  New  Mexico. 

The  attitude  measuring  system  of  this  rocket  consisted  o.f  a set 
of  triaxlal  Schonstedt  600  milligauss  magnetometers  and  a Bayshorc 
System  Lunar  Sensor  Model  number  LS-ll-DR-2.  Two  magnetometers  which 
were  perpendicular  to  the  rocket  axis  were  in  the  same  plane  as  the  lunar 
sensor  but  45°  out  of  phase  (refer  to  Figure  1)  therefore,  the  readings 
of  magnetometers  had  to  be  phase  shifted  so  that  attitude  system  would 
give  simultaneous  readings  of  Earth's  magnetic  field  and  the  moon. 

The  lunar  sensor  angular  output  was  saturated  for  nearly  the  entire 
flight  and  only  functioned  for  a few  discrete  points  on  the  descent 
portion  of  the  flight,  but  the  times  of  the  lunar  readouts  were  determined 
accurately.  To  calculate  the  angle  between  the  rocket  axis  and  the 
lunar  vector,  techniques  not  included  in  a normal  reduction  of  rocket 
attitude  data  were  used  with  acceptable  results.  The  approach  undertaken 
incorporated  modifications  of  existing  methods  with  additional  criterion 
studies  [2]  [3] . The  primary  data  needed  was  reliable  magnetometer  data. 
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both  longitudal  and  lateral  from  which  the  approximate  angular  velocity 
of  precession  and  spin  rate  could  be  determined,  and  available  lunar 
information.  The  above  approach  assumed  the  rocket  to  be  well-behaved 
for  some  portion  of  the  flight. 


The  basic  coordinate  system  used  in  this  report  is  a true  north, 


e , e 
6 r 
c c c 


east,  and  vertical  cartesian  coordinate  system  where  e . 
it  a system  of  orthonormal  vectors  with  tangent  to  the  meridan  circle 


c 

through  the  launching  tower  and  pointing  to  true  north,  e tangent  to 

^c 

the  circle  of  latitude  through  the  launching  tower  and  pointing  east  of 

north , and  S be  determined  by 
* r . 


This  system  will  be  referred  to  in  this  report  as  the  north,  cast,  and 
vertical  system  (refer  to  Figure  2). 


i 

[ 

! 

! 

> 

t' 

\ 


- 2 - 


SECTION  1 


ATTITUDE  DETERMINATION  AND  REFINEMENTS  FROM  GYROSCOPIC  DATA 

1 . 1 MAR*  AND  MIDAS  GYROSCOPIC  ATTITUDE  SYSTEMS 

The  Whittaker  Corporation's  Miniature  Attitude  Reference  System  (MARS) 
and  Space  Vector  Corporation's  Miniature  Inertial  Digital  Attitude  System 
(MIDAS)  both  provide  gyroscopic  outputs  for  the  determination  of  rocket 
attitude.  Both  systems  reference  roll,  pitch  and  yaw  motion  in  the  same 
manner;  i.e.,  the  roll  axis  coincides  with  ine  rocket  axis  at  launch 
and  the  orientation  of  a gyro  reference  notch  determines  the  location  of 
the  yaw  axis.  The  pitch  axis  is  fixed  as  the  cross  product  of  the  yaw 
axis  with  the  roll  axis. 

Although  the  gyro  notch  is  normally  set  in  the  vertical  plane 
containing  the  rocket  axis  at  launch,  it  may  be  rotated  to  any  orientation 
perpendicular  to  the  rocket  axis.  Since  the  accuracy  of  the  yaw  and 
p tch  depends  upon  the  orientation  of  the  notch  [1] , it  is  essential  that 
its  exact  position  be  known. 

Both  the  MARS  and  MIDAS  systems  employ  a roll -stabilized  platform, 
using  two  two-degrees-of-freedom  gyroscopes.  The  gyros  are  supported 
by  gimbals  arranged  on  the  platform  so  that  the  gimbals  are  free  to 
rotate  about  the  roll,  pitch  and  yaw  axes.  Initially,  the  orientation 
of  the  platform  is  fixed  at  launch.  When  the  vehicular  motion  begins  to 
change  the  orientation  of  the  platform,  precession  of  the  gyroscopes 
is  detected  by  pickoffs  on  the  corresponding  output  axes.  Signals 
transmitted  to  a SERVO-mechanism  then  activate  a motor  which  restores 
the  platform  to  its  initial  attitude.  Thus,  the  stabilized  platform 
retains  its  orientation  in  inertial  space,  providing  a fixed  coordinate 
reference  system  for  roll,  pitch  and  yaw  measurements. 

The  outputs  from  the  MARS  and  MIDAS  system  differ  considerably,  and 
consequently  different  techniques  are  used  to  recover  angular  roll, 
pitch  and  yaw.  The  MARS  system  employs  potentiometers  to  convert  roll, 
pitch  and  yaw  displacements  to  analog  voltages,  having  a range  of 
approximately  0-5  volts.  Each  5-volt  span  represents  some  degree 
segment  having  a predetermined  length.  A sudden  change  in  the  voltage 
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(from  5 to  0 volts  or  from  0 to  5 volts)  is  called  a swap,  and  it  is 
only  through  the  location  of  swaps  that  the  correct  degree  span  can  be 
determined  from  the  MARS  data  [i]. 

l;or  a vehicle  equipped  with  the  MIDAS  system,  the  roll,  pitch 
and  yaw  measurements  are  digitized  by  means  of  optical  encoders.  As 
a result  the  one-to-one  correspondence  between  angular  displacements 
and  digital  codings  satisfies 


where 

(j)  Is  the  angular  displacement  in  degrees, 

n is  the  digital  coding 

and  1024  counts  represent  360  degrees.  No  bias  voltages  ai’e  needed 
for  conversion  as  with  the  MARS  system  since  the  MIDAS  encoders  auto- 
matically represent  the  displacement  from  the  uncaged  position  of  the 
gyro  which  fixes  the  coordinate  reference  system. 

1 • “ Attitude  Control  System 

For  i he  Space  Vector  Attitude  Control  System  (ACS),  the  orientation 
of  the  gy' o notch  is  not  used  as  a r:ferencc.  Instead,  the  orientation 
of  the  "Top  Dead  Center"  (TDC)  supercedes  all  other  references,  and 
all  side  probe  measurements  are  based  upon  the  TDC  position. 

1.2.1  ACS  Coordinate  Reference 

The  ACS  coordinate  reference  is  displayed  in  Figure  3.  The  yaw, 
pitch  and  roll  references  duplicate  those  for  the  MARS  and  MIDAS 
systems  except  that  gyro  notch  position  is  not  available.  The  TDC 
reference  as  used  for  the  ACS  system  is  displayed  in  Pigure  4. 

1.2.2  ACS  Output  Measurements 

Three  models  for  the  ACS  were  encountered  in  attitude  determination  - 
II 10370,  // 1 0380 , II 10390.  However,  for  data  reduction,  models  10370 
and  10380  were  treated  alike.  Model  10370  (10380)  output  roll 
voltage,  pitch  course  and  fine  voltage,  yaw  course  and  fine  voltage. 
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Model  10390  output  only  roll  voltage-pitch  coarse  and  yaw  coarse  voltage. 
The  conversion  to  angular  measurements  was  made  with  calibration  curves 
which  provided  the  magnitude  of  the  angle. 

The  following  Table  1 summarizes  the  ranges  of  the  typical  calibration 
curves  for  the  various  models.  It  should  be  noted  however,  that  linear 
interpolation  was  not  possible  with  all  conversions. 


MODEL 


10370 

10380 

10390 

VOLTS 

DEGREES 

VOLTS 

DEGREES 

VOLTS 

DEGREES 

0 

Pitch  Coarse 

(.54,4.36) 

0 + 80 

(.53,4.33) 

0 + 80 

(-.05,4.64) 

0+80 

u 

Yaw  Coarse 

(1.37,3.56) 

0 + 45 

(1.37,3.56) 

0 + 45 

(.95,  4.14) 

0+50 

T 

P 

Roll 

(0,  4.63) 

0 + 170 

(0,  4.62) 

0 + 170 

(0,  5) 

0+45 

U 

Pitch  Fine 

(-.54,4.63) 

0+8 

(-.53,4.78) 

0 i 7 

- 

- 

T 

Yaw  Fine 

(-.52,4.94) 

0 + 7 

(-.51,4.85) 

0 + 7 

- 

" 

TABLE  1 


The  sign  reference  for  outputs  can  be  determined  from  Figure  3 along 
with  the  following  Table  2. 


MODEL 


10370  5 10380  j 

10390 

ROLL  CIV 

PITCH  DOWN 

YAW  RIGHT 

ROLL  CW 

PITCH  DOWN 

YAW  RIGHT 

Pitch  Coarse 

- 

decreasing 

- 

- 

decreasing 

- 

Yaw  Coarse 

- 

- 

increasing 

- 

- 

increasing 

Roll 

increasing 

- 

- 

decreasing 

- 

- 

Pitch  Fine 

- 

increasing 

- 

- 

- 

- 

Yaw  Fine 

“ 

decreasing 

- 

- 

- 

TABLE  2 
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1*3  Analysis  and  Refinement  Techniques 

The  development  of  analyses  to  calculate  the  attitude  of  a rocket 
from  gyroscopic  roll,  pitch  and  yaw  data  are  discussed  in  the  following 
subsections.  Since  telemetered  data  may  possess  inherent  noise,  data  filtering  and 
smoothing  may  be  required  to  provide  continuous  output.  To  substantiate 
the  curve  fitted  data  in  the  smoothing  process,  statistical  studies  are  ; 

made  to  compare  measured  and  predicted  data  discrepancies.  ( 

*•3.1  Analysis  for  Determination  of  Rocket  Attitude  ; 

*.ccording  to  [1] , the  longitudinal  axis  of  the  vehicle  can  be 
expressed  as  the  vector  e": 

r 1 

A /v  a 

= X cosy  cosp  + Y siny  + Z cosy  sinp  : 

where  y and  p are  the  yaw  and  pitch  respectively  and  the  body  axes  X, 

Y,  Z are  in  the  direction  of  the  gyro  roll,  pitch  and  yaw  axes. 

Given  the  elevation  0 and  azimuth  <j>  of  e"  in  the  local  North 

r 

East  and  Vertical  system  SQ  , e § at  launch,  the  unit  vector 
* * * c c 

X,  Y,  Z can  be  expressed  as  linear  combinations  of  the  earth-based- 

system  eg  , 6.  , e so  that  we  have 
c ?c  rc 


where  the  a^  terms  (i,j  = 1,2,3)  are  the  initial  direction  coefficients 
of  the  gyro  axes  at  launch. 

According  to  equation  (5)  of  reference  [1],  we  can  define  a system 
of  orthonormal  vectors  e"  c",  e"  at  any  time  in  flight  as  a linear 

A A 1 A X ilm 

combination  of  the  X,  Y,  Z axes.  The  elevation  0 and  azimuth  <j>  of  the 
rocket  axis  or  any  vector  perpendicular  to  the  rocket  axis  can  now  bo 
determined  by  (9)  and  (10)  described  in  reference  [1]. 


/ 
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For  side  probe  data,  the  orientation  of  the  probe  must  be  defined 
in  the  e^,  e^,  e^'  system;  i.e.,  for  P in  the  direction  of  the  probe 

A 

P = e"  cosX  + e”  cosy  sinX  + elj  siny  sinX 

/\ 

where  X is  the  angle  between  e"  and  P,  y is  the  angle  between  e'1 

a r i 

and  the  projection  of  P in  the  plane  of  e^  and  e'^. 

Since  e",  e"  e£  can  already  be  expressed' in  the  local  Sg  , 
r c 

e,  , c system,  only  a rotation  matrix  involving  X and  y is  necessary 
^c  rc 

to  express  P as 

/N 

P = ert  cos0  cos4»  + S,  cos0  sind>  + 6 sin0  . 

9c  P P 4»c  P > rc  P 

The  angle  0^,  measured  positive  up,  is  the  elevation  of  the  probe 

axis  above  the  horizontal  plane  of  §-  and  e.  , and  <f>  is  the  azimuth 

°c  ^c  ** 

of  the  probe  axis  measured  from  true  North  positive  east.  Throughout 

the  remainder  of  this  report,  all  such  angles  0 and  <j>  describing  the 

orientation  of  a probe  or  sensor  axis  in  the  local  , §.  , $ 

c ^c  c 

system  will  reference  elevation  and  azimuth  respectively. 

A 

Once  the  components  of  P in  the  above  system  are  determined  and 
the  components  of  the  unit  vector  in  the  direction  of  the  rocket  flight 
are  given  by 

v * 6e  vj  + % v2*  % v3 

c rc  c 

/S  A 

then  the  attack  angle  i p between  P and  V is  found  from 


i p = cos"1  (VjCosOp  cos<J>p  + V2Cos0pSimJ>p  + v^sindp) 


Refer  to  Figures  (5  - 11)  for  typical  angles  of  attack. 
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1.3.2  Refinement  Techniques  for  Final  Outputs 

Due  to  noise  or  in-flight  calibrations  often  times  present  in  the 
recorded  data,  filtering  procedures  are  modified  or  developed  as  required 
to  provide  continuous  final  output  for  the  attitude  determination 
system.  The  primary  routines  used  are  Fourier  Series  and  nth  degree 
polynomial  approximations. 

Typical  examples  of  the  above  recorded  data  problems  are  encountered  in 
the  conversion  of  the  ACS  coarse  (and  fine  when  available)  yaw  and  pitch 
data  (see  Section  1.2.2)  to  angular  measurements.  The  random  noise 
discontinuities  in  these  converted  outputs  cause  the  attitude  determination 
results  to  be  unusable  in  the  analysis  of  probe  data  requiring  continuity 
of  vehicle  motion  (refer  to  Figures  12  and  13). 

1 .3.2.1  Fourier  Series  Approximation 

To  provide  smooth,  continuous  and  acceptable  attitude  information 
during  the  well-behaved  areas  of  a particular  vehicular  flight,  the 
coarse  (and/or  fine  when  applicable)  yaw  and  pitch  data  can  be  predicted 
by  the  Fourier  Series 

Y = + s (acosS-  + bsin—)  (1.3-1) 

, II  Li  n L# 

n=l 

This  approach  is  used  due  to  the  periodic  nature  of  the  data. 

For  most  rocket  flights,  quick  convergence  for  (1.3-1)  is  obtained 
with  the  approximation 

Y = A + Bt  + C sinwt  + D coscot.  . (1.3-2) 

The  addition  of  the  linear  term  to  the  1st  three  terms  of  (1.3-1)  is 
needed  to  account  for  an  occasional  linear  shift  in  the  converted 
data.  A preliminary  rotation  rate  or  angular  velocity  w is  selected 
from  a study  of  the  oscillograms  of  the  raw  data,  and  is  further  refined 
by  an  option  within  the  fitting  routine. 

1.3. 2. 2 Nth  Degree  Polynomial  Approximation 

During  regions  in  which  a vehicle  was  not  well-behaved,  yaw  and 
pitch  outputs  are  separated  into  discrete  time  intervals.  These 
intervals  are  then  curve  fitted  with  polynomials  up  to  the  20th  degree 
when  necessary.  The  nth  degree  polynomial  routine  selects  the  best 
fit  for  the  data  calculating  an  RMS  value  (see  Section  1.3.3)  between 
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measured  data  and  data  calculated  from  eaJi  of  the  polynomial  approximations. 
The  minimum  RMS  value  determines  the  best  fit  in  the  specified  interval. 

These  techniques  were  applied  with  acceptable  results  to  A17. 110-1  during 
regions  of  ascent  and  descent  with  the  vehicle  was  not  well-behaved. 

1.3.3  Statistical  Studies 

For  each  of  the  data  refinement  techniques  discussed  in  Section  1.3.2, 
RMS  values  are  computed  for  the  polynomial  or  series  approximation  to 
the  measured  output.  However,  a small  RMS  value  can  sometimes  be  mis- 
leading; i.e.,  slight  unwarranted  wave  patterns  may  be  created  by  (1.3-2) 
in  the  predicted  data  whereas  plots  of  the  converted  yaw  and  pitch  data 
indicate  a strong  possibility  of  a linear  fit  to  the  data.  In  cases  like 
described  in  Figure  14c,  additional  studies  of  the  correlation  coefficient 
p are  made  to  check  whether  a linear  fit  or  equation  of  the  form  (1.3-2) 
would  best  represent  the  data.  Computation  of  the  RMS  value  together 
with  the  correlation  coefficient,  when  required,  form  an  .integral  part  of 
the  reliability  assessed  to  the  fitting  routines  described  in  Section  1.3.2. 

1.3. 3.1  RMS  Value  with  Application 

The  RMS  (root  mean  square)  deviation  of  a set  of  N values  is  defined 
to  be  the  square  root  of  the  mean  of  their  squares;  that  is 

-&  (X.-X)2 

RMS  (1.3-3) 


where  is  the  difference  between  a measured  and  time  associated  predicted 
data  value,  and 


X 


To  simplify  computation,  (1.3-3)  can  be  rewritten  as 


RMS 


(1.3-4) 


The  RMS  deviation  can  also  be  referred  to  as  the  standard  deviation  0. 
For  a normal  distribution,  a +_  3a  deviation  accounts  for  =99.7%  of  the 
measured  data.  Therefore,  tolerance  levels  are  set  on  input  data  to 
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ensure  that  all  predicted  data  be  within  the  specified  a deviation 
utilising  the  results  of  (1.3-4). 

For  the  vehicle  A1 8. 006-4,  (1.3-2)  was  used  to  generate  yaw  and 

pitch  data  from  80  to  375  seconds  after  launch.  This  function  took  into 

account  the  sinusoidal  nature  of  the  output  as  well  as  the  slight  linear 

displacement  of  the  data.  Using  a sample  rate  of  20  points  per  second 

with  an  acceptance  level  of  1°  between  data  points  as  input  to  the  software 

2ir 

routine,  output  showed  the  unknowns  in  (1.3-2)  to  be  w = yj  7273 
the  pitch  coefficients 

A = 6.863463  B = -.0026176  C = .132396  D = .090675 

and  the  yaw  coefficients 

A = .6111509  B = .0020662  C = -.030052  D = .179948. 

Generating  data  with  the  above  coefficient  values,  the  RMS  values  for 
pitch  and  yaw  were  shown  to  be  0.565°  and  0.288°  respectively. 

1.3.3. 2 Correlation  Coefficient  with  Application 

For  data  points  (t^,  y^),  the  correlation  coefficient  p is  defined 
as 


N 

N N 

( S t±)  (2  y.) 

1=1 

i=l  i=l 

- N N 

N N 

NS  t . 2 - (2  t.)2 
i=l  i=l 

N 2 y?  - (2  y-)2 
_ i=l  i=l 

) 

p may  take  on  any  value  between  -1  and  +1  inclusive;  and  p close  to 
-1  implies  a strong  negative  correlation  whereas  p close  to  +1  implies 
a strong  positive  correlation,  p close  to  0 implies  a weak  correlation. 
Several  examples  of  the  strengths  of  linear  correlation  [5]  are  given  in 
Figure  14.. 

When  displaying  the  pitch  fine  and  yaw  fine  outputs  for  A 18. 11 6-1, 
a strong  negative  correlation  for  each  was  found  from  100  to  350  seconds 
after  launch.  Application  of  (1.3-2)  to  converted  pitch  fine  and  yaw 
fine  data  showed  C - 0 and  D = 0 as  were  expected.  The  expression  for 
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tlie  predicted  pitch  data  in  degrees  was 


PITCH  = - . 00015545T  + 7.021074,  100  _<  T < 350 

with  a resulting  RMS  of  0.06°.  This  RMS  value  was  sufficiently  low 
to  allow  neglecting  other  possible  vehicle  motions.  Similarly,  the 
expression  for  the  predicted  yaw  data  in  degrees  was 

YAW  = .000217S66T  + .4219051,  100  _<  T < 350 

with  a resulting  RMS  of  0.05°.  Again,  the  RMS  value  was  well  within 
acceptable  error  limits  so  as  to  allow  neglecting  other  possible  motions. 

1 .4  Vehicle  Studies 

To  corroborate  the  initial  orientations  of  on-board  sensors, 
measured  and  attitude  predicted  magnetometer  pitch  angles  are  compared 
[1].  A software  package  was  developed  to  compare  the  phase  relationship 
during  well-behaved  areas  of  a vehicle's  flight. 

This  package  involves  a least  squares  study  on  the  magnetometer 

data  (refer  to  Appendix  A for  detailed  description).  The  output 

produced  includes  a table  of  maximum  (minimum)  times  T measured 

data  values  in  degrees  and  associated  AT  values  for  the  corrected 

meas 

magnetometer  outputs.  Also  included  are  a table  of  maximum  (minimum) 

times  T , predicted  data  values  in  degrees  and  associated  AT 

gen  gen 

values  for  the  predicted  magnetometer  output  from  the  attitude  data. 

The  two-time  dependent  tables  are  merged  for  all  consecutive  time 
values  of  each  array  and  the  phase  shift  between  the  associated  values 
are  printed. 

For  T < T + AT  , the  phase  shift  x between  the  raw 
gen  — meas  meas  A 

and  the  predicted  magnetometer  data  is 


X 


T 

gen 


- T 

meas 

7\T 


meas 


x 360. 


l:or  0 <_  x £ 180  the  predicted  data  is  lagging  or  coincident  with  the 
measured  magnetometer  data.  For  180°  < x < 360°,  the  measured  data  is 
said  to  lag  the  predicted  by  360°  - x*  The  following  Figure  is  depicts 
the  relationships  just  discussed: 


PHASE  RELATIONSHIPS 


l-igure  15 


fn  this  figure,  the  predicted  data  from  the  attitude  is  lagging  the  raw 


data  since  we  are  assuming  that  (T 


- T ) < L AT 
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SECTION  2 


ATTITUDE  SYSTEM  FOR  A09. 209-1 


The  primary  attitude  sensing  devices  on  board  consisted  of  a 
triaxial  set  of  magnetometers  and  a lunar  sensor.  Also  aboard  was 
a set  of  photomoters  which  provided  approximate  lunar  information. 
From  these  devices  (if  they  function  properly),  it  i;'  possible  to 
calculate  the  attitude  for  any  probe  at  any  time  during  the  flight. 


f 


2 . 1 Lunar  Sensor 

The  Bayshore  lunar  sensor  had  its  sensitivity  set  to  read  out 
a wide  range  of  high  lunar  phase  angles  with  a resolution  of  2°  and 
an  accuracy  of  +1°.  The  sensor  was  mounted  in  such  a way  that  it 
reads  out  only  when  the  rocket  axis,  the  moon,  and  the  sensor  axis  are 
all  in  the  same  plane  (once  per  spin  of  the  rocket).  The  sensor  output 
which  is  in  the  form  of  a digital  code  can  be  converted  directly  to 
lunar  aspect  angles  (refer  to  Figure  16).  The  angles  this  sensor  was 
capable  of  reading  ranged  from  +70°  to  -26°  as  measured  from  the  center 
line  of  the  sensor. 


FIELD  OF  V1EIV  OF  LUNAR  SENSOR 


T 

+Rockct 

Axis 


Figure  17 
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2.2  Photometers 


On  hoard  this  vehicle  was  a set  of  three  solar  blind  photo- 
multipliers with  UV  interference  filters  with  a limited  field  of  view. 

This  set.  of  photometers  was  mounted  in  such  a manner  that  the  projections 
of  the  three  photometers  on  the  plane  perpendicular  to  the  rocket 
axis  were  coincident  with  projection  of  the  lunar  sensor,  L.e.,  the 
azimuth  angles  in  the  horizontal  plane  w;ro  equal.  However,  only  one 
of  the  photometers  was  aligned  with  the  lunar  sensor  and  this  photometer 
failed  to  function  properly.  The  others,  one  with  a 2750  Angstrom  filter 
which  was  mounted  up  5°  from  the  horizontal  plane  and  the  other  with 
2600  Angstrom  filter  which  was  mounted  up  10°  from  the  horizontal  plane 
functioned  properly,  (see  Figure  1). 

The  output  from  these  two  photometers  provided  some  important 
lunar  information  since  they  indicate  maximum  intensities  when  they  sense 
the  moon. 

2.5  Magnetometers 

The  magnetometer  system  on  board  was  a triaxial  flux-gate  magnetometer, 
specifically  designed  for  use  in  rocket  aspect  measurement,  systems . 

A field  component,  within  the  range  of  + 600  inilligauss  is  converted 
to  an  analog  voltage  defined  by  the  equation  of  the  type 

E = 2.40  + .004  H cos<|> 

where  1:  is  the  output  in  volts,  II  is  the  ambient  field  in  milldgauss 
and  <{»  is  the  angle  between  the  magnetic  field  vector  and  sensor’s 
positive  magnetic  axis.  Figures  18  and  19  show  the  specifications 
and  the  vendors  calibrations  for  this  system.  This  calibration  is 
done  on  the  instrument  not  under  simulated  flight  conditions. 

2.4  Errors  in  Magnetometer  Readouts 

Since  these  magnetometers  sense  the  total  ambient  field,  the 
need  of  reliable  magnetometer  information  requires  that  all  permanent, 
stray,  and  induced  fields  should  be  determined  as  accurately  as  possible. 
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On  this  particular  flight,  there  was  a Siemens  Brushless  DC  motor 
with  four  permanent  magnets  which  spun  at  approximately  3,000  RPM. 

This  motor  was  located  30°  off  the  longitudinal  axis  and  in  a plane 
parallel  to  door  1,  (refer  to  Figure  1).  This  motor  modulated 
the  voltage  output  of  the  three  magnetometers.  It  created  in  addition 
to  normal  permanent  and  induced  fields,  two  effects;  first,  it  created 
a permanent  magnetic  field  on  all  three  axes  which  were  modulated  at 
the  spin  frequency  of  the  motor.  This  effect  on  each  axis  is  clearly 
demonstrated  in  Figures  (20-22)  where 

X axis  « .112  volts  - 26.88  milligauss 

Y axis  - .147  volts  - 35.28  milligauss 

Z axis  - .3  volts  = 72  milligauss 

Secondly,  since  the  axis  of  the  motor  was  not  aligned  parallel  to  a 
magnetometer  axis,  it,  too,  was  spinning  at  the  vehicle  spin  rate. 

This  fact  coupled  with  the  absence  of  preflight  calibration  of  the 
magnetometers  under  simulated  flight  conditions,  makes  orthogonality  in 
the  measurements  of  the  magnetometers  almost  impossible. 

2.5  Permanent  and  Induced  Magnetic  Fields 

Since  there  was  no  preflight  calibration  of  this  set  of  magnetometers 
under  flight  conditions,  we  were  unable  to  determine  all  the  induced 
magnetic  fields  on  each  magnetometer  axis.  Instead,  the  following 
approach  was  utilized. 

To  minimize  magnetometer  output  error,  previously  written  routines* 
have  been  implemented  to  correct  for  bias  changes  and  induced  field 
contributions  in  a mutually  perpendicular  three  magnetometer  system. 

Given  the  bi  cor.  ribution  parameters  c^,  i=l,2,3  for  the  X,  Y and 
Z magnetometer  outputs,  respectively  and  the  theoretical  field  H,  we 
can  form  the  equation 


H2  = A,[X  + a.]2  + A-[Y  + a-]2  + A,[Z  + aj2. 
I1  mg  1J  2l  mg  2J  3l  mg  3J 

X , Y , Z are  the  results  of  vendor  supplied  calibration  curves 
mg’  mg’  mg 


*LS0.MAG,  Paul  Pruneau,  Space  Data  Analysis  Laboratory,  Boston  College 
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which  convert  raw  magnetometer  output  to  milligauss.  The  least  squares 
induced  field  constants,  i=l>2,3  are  determined  by  forming  the 

sum 


I 

n 


r[H2  - A,  (X  +CC.)2  - 
* 1 1 mg  V 


A~(Y  +ot_)' 
2^  mg  2J 


A-  (Z  +a..)“J ' 
3V  mg  3 J 


and  setting 


0 


i= 1,2,3 


This  will  form  a system  of  three  equations  with  three  unknowns  and 
a unique  solution  can  be  found  for  the  A^  constants  so  long  as  the 
determinant  of  the  coefficients  docs  not  vanish. 

The  software  package  iterates  upon  the  ol  values  for  a given 
iteration  interval  and  determines  the  A.  constants  for  each  iteration. 

l 

An  RMS  (Root  Mean  Square)  is  then  computed  between  the  measured  total 
magnetic  field  and  the  theoretical  total  field.  Those  values  of  cu 
and  associated  +V\7  values,  i=l,2,3  for  which  the  RMS  is  a minimum  are 
the  permanent  and  induced  field  corrections  to  the  magnetometer  output. 
A flow  chart  of  the  logic  contained  in  this  routine  as  displayed  in 
Figure  23.  This  method  gives  the  best  average  inducing  effect  on  each 
axis,  but  it  docs  not  give  a true  representation  of  the  induced  field.;. 

A true  representation  of  the  outputs,  where  the  actual  magnetic 
field  on  each  magnetometer  axis  was  known,  would  be  given  by  the 
system 


where 


X'  , Y*,  Z» 
m’  m’  m 


V Yt>  zt 


x'  - 

+ 

b.  X. 

+ b-,Y 
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+ X 
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are  the  actual  magnetometer  outputs  c'>- 1 ilnin  , 
all  types  of  magnetic  influence 

are  the  true  magnetic  fields  that  the  magnetometers 
should  be  reading 


X Y Z 
perm’  perm’  perm 


b.  - b 
J 9 


are  the  permanent  fields  on  each  axis 
are  the  induced  fields 
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Since  accurate  permanent  fields  (bias)  were  needed  for  later  analysis, 
additional  studies  were  performed  as  outlined  below. 

1.5.1  Data  Sample  Rate  and  Interval  Section 

In  order  to  accurately  determine  the  times  of  maximum  magnetometer 
output  (see  Appen  lix  A),  smooth  magnetometer  readings  are  required. 

Normal  data  transmission  does  not  provide  such  smooth  output,  so  to 
allow  for  proper  usage  of  the  routines  described  in  Appendix  A,  a 
curve  fit  was  put  through  the  bias  values  of  the  magnetometer  data. 

The  samples  selected  for  this  fit  were  tak.;n  from  the  well-behaved 
portion  of  flight  - i.e.,  from  approximately  40  kms.  on  ascent  to 
approximately  40  kms.  on  descent.  Since  the  rotation  rate  of  the 
magnetometer  output  was  - .175  seconds  and  every  other  point  of  the 
available  333  samples  per  second  for  each  magnetometer  was  used,  we  had 
available  166  points  per  second  cr  29  points  per  rotation. 

Three  120  second  time  intervals  [65,185),  [185,305),  and 
[305,425)  were  used  to  obtain  a separate  bias  approximation  B for  each 
interval  from  the  equation 

Z(VrB)  = 0 

where  V.  represents  the  data  sample.  This  equation  reduces  to 

SVi 

B=  ~~rT~  ' 

In  the  120  second  time  interval,  a simple  calculation  shows  that  the 
number  of  data  points  n is  20,000. 

2.5.2  Polynomial  Fit 

The  standard  bias  equation  as  a function  of  time  is  approximated 
by  a first  degree  polynomial.  However  to  obtain  greater  accuracy 
with  higher  order  terms  in  the  fitting  routine  results.,  the  value  of  the 
bias  B(t)  as  a function  of  time  t was  equated  with  a 2nd  degree  polynomial 
to  give 

B (t)  = at2  + bt  + c.  vv 

^ v ; 

However,  from  the  discussion  above  wc  can  approximate  B(t)  by 


in  each  of  the  intervals  studied.  This,  in  turn,  leads  to  the  system  of 
equations 


at*  + bt^  + c 
at2  + bt2  + c 
at|  + btg  + c 


where  t . , j = 1 ,2 ,3  is  taken  as  the  mean  time  in  the  jth  interval  and 


•V* 

i=l 


3 = 1, 2, 3 


As  tiie  value  of  n increases,  the  maximum  error  that  occurs 
in  the  approximation  of  15^  decreases.  This  fact  can  be  verified  by 

SVi 

c Iculations  on  the  expression — — . 


This  procedure  to  determine  the  bias  equation  as  a function 
of  time,  forms  part  of  the  procedure  used  on  all  vehicle  flights  when 
examining  phase  shifts  between  measured  and  predicted  data. 

2 . 6 Magnetometer  Phase  Shifting 

To  properly  utilize  the  attitude  determination  technique  mentioned 
in  [2][3]it  is  necessary  to  have  simultaneous  readings  of  the  Harth's 
magnetic  field  and  the  moon  from  attitude  sensing  devices  with  the  same 
orientation  on  the  vehicle.  On  this  particular  vehicle,  both  the  lunar 
sensor  and  both  the  lateral  magnetometers  were  mounted  normal  to  the 
axis  of  the  rocket,  but  they  were  out  of  phase  by  45°  (refer  to  Figure 
1). 

To  properly  align  a magnetometer  axis  parallel  to  the  lunar 

/■>  A 

sensor  axis,  we  begin  by  defining  the  unit  vectors  Y and  Z along  the 
Y and  Z magnetometers  axes  respectively  as 


Y = cosO  cos*  * 6,  cosO  sind)  +6  sinO 

0C  y V <PC  y y rc  y 

A 

Z = 6.  cos0  cosd>  + 6,  cos0  sin<j>  + 6 sinO  . 

0 z Tz  <f>  z rz  t z 

c rc  c 
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Then  the  unit  vector,  U,  along  the  lunar  sensor  axis  and  in  the  plane 
of  Y and  Z,  can  be  expressed  as 


/S  A 


U = Y cosx  + Z sinx 

where  x is  the  angle  between  the  lunar  sensor  axis  and  Y magnetometer 
axis . 

Then  the  scalar  product  of  the  Eartn's  unit  magnetic  field  vector  H 
with  tl  where 


H = e.  x + c , y +6  z 
0 m <b  /m  r m 
c Yc  c 


is  given  by 


cosy  = (H*U)  = (H»Y)cosy  + (H*Z)sinx 


where  (H*Y)  and  (H«Z)  are  equal  to  the  actual  normalized  magnetometer 

outputs  H and  H . 

y z 


Then 


Hu  = cosy  - ii^cosy  + 1{zsinx 

is  the  expression  that  represents  the  cosine  of  the  angle  between  the  magnetic 
field  and  the  lunar  sensor  axis. 

Accurate  knowledge  of  the  bias  of  this  generated  magnetometer 
data  (or  the  times  when  this  magnetometer  is  reading  approximately 
zero  magnetic  field)  is  helpful  in  the  attitude  analysis.  Therefore, 
the  numerical  methods  previously  discussed  in  sections  (2.5.1)  and 
(2.5.2)  were  applied  to  this  data. 
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SECTION  3 


ATTITUDE  ANALYSIS  FOR  THE  WELL-BEHAVED  PORTION  OF  THE  FLIGHT 


The  basic  attitude  determination  analysis  for  this  vehicle  uses 
a combination  of  available  lunar  information  and  corrected  magnetometer 
data.  The  available  lunar  information  was  obtained  through  interpreta- 
tion of  the  behaviorial  pattern  of  the  photometer  data  at  selected  times, 
through  the  accurate  times  of  the  lunar  sensor  readout  (even  though 
the  sensor  angular  output  was  saturated  for  the  most  of  the  flight),  and 
through  techniques  which  allow  the  lunar  angle  to  be  determined  when 
specified  criteria  are  fulfulled. 

3. 1 Determination  of  the  Elevation  and  Azimuth  of  the  Rocket  Axis  at  a 
Lunar  Readout  Time 

AAA  A 

Let  X,  U,  W be  a orthonormal  system  with  X,  a unit  vector  in 
direction  of  the  rocket  axis,  U,  a unit  vector  perpendicular  to  the  rocket 

AAA 

axis  at  the  lunar  sensor,  and  X x U = W.  Since  we  know  the  lunar  sensor 

A A 

reads  out  only  when  the  moon  is  in  the  plane  of  X and  U, 

M = M X + M U (3.1-1) 

x u 

and  representing  the  magnetic  field  in  this  vehicle  system 

II  = H X + II  U + II  W.  (3.1-2) 

X u w 

In  (3.1-1)  and  (3.1-2)  as  well  as  in  the  following  discussions,  M and  II  are 
assumed  to  be  unit  vectors.  The  scalar  product  of  the  vectors  II  and  M is  then 
given  by 

(II»M).  = II  M + M H . (3.1-3) 

v '1  X X u u 

The  scalar  product  can  also  be  derived  from  empiracal  sources  where 


and 


M = m^  + 
c 1 


e<j>  m2  + er  m3 
rc  c 


(3.1-4) 


H = 6.  x + §,  y + c z . 
0„  m <\>  r m 

c rc  c 


(3.1-5) 
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Then 


(h-M)2  = m,xB  * »2ym  * .J*, 


m 


(3.1-6) 


From  (3.1-3)  and  (3.1-6) 


(H-M)2  « HxMx  + MuHu 


(3.1-7) 


and  if  the  magnetometer  parallel  to  the  axis  reads  a zero  magnetic  field 
(refer  to  section  2.6)  at  the  time  of  a lunar  sensor  readout  then  from 
(3.1-7),  it  is  possible  to  determine  an  angular  output  from  the  lunar 
sensor  since 

(H«M)0 


and 


Mx  = — 


Mu  = sin(cos-1 (Mx)) 


The  elevation  0 and  the  azimuth  <b  , of  the  rocket  axis  from  the 
x x’ 

equations  on  page  8 of  [3] . 

(Hx-(H‘M)Mx)zm  + (Mx-(H-M)Hx)m3  + (HMX)  CyWnP 


sin0 


1 - (H*M): 


(3.1-8) 


(U  - (H*M)M  )y  + (M  - (H»M)H  )m-  + (HMX) (x  m -z  m ) 

. x xm  x a l m o in  i i o') 

n<J>x  = (Hx-(H*M)Mx)xm  + (Mx-(H*M)Hx)mi  + (1WX)  (z^-m^f 


where  from  (3.1-1)  and  (3.1-2) 

✓v 

(HMX)  = -HwMu. 


3.2  Determination  of  the  Axis  of  Precession 

After  an  orientation  of  the  rocket  axis  during  the  well-behaved 
portion  of  the  flight  has  been  determined,  the  equations  in  [2] 
page  15,  were  applied  to  determine  the  axis  of  precession. 

Then  0^,  the  elevation  of  the  axis  of  precession,  and  <f>^,  the 
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azimuth  are  given  by 


0p  = sin 


■i 


| (cosgH-cosY1cosa)zm+Ccosa-cosY1cos3H)sinex+(HeNer) (HeN§r  ) 


sxnY, 


(3.2-] 


tan"1 


(cos3H-cosYjCOsa)ym+  (cosa-cosY1cos3^)cos0xsind>x+  (He^e  ) (He^e , ) 

c 

(cosp^-cosYjCosoOx^  (cosa-cosY^os3H)cos0xcos<bx+(HeNer)  (He~£~  )" ' 


(3.2-: 


where 


an  orientation  of  the  rocket  axis  X is  given  by 

ert  cosO  cosd)  + B,  cos0  sincb  + $ sin0 
0 x Yx  d)  x Yx  r x 
c rc  c 


e^  = the  axis  of  precession 

= en  cos0  cosd)  + B,  cos0  sin<j>  + 3 sin0 
6C  P P 4>c  P P rc  p 


cosYj  = (H*eN) 

cosa  = % cone  angle  = (e^e^) 

cos3H  = the  angle  the  magnetic  field  vector  makes  with 
axis  of  precession  = (H*e^) 

(IIcNer)  =i  ^si n2"a  sin23f^-  (cosYj  -cosS^cosa) 2 

(He.,§  ) = y cos0  cosd>  -x  cos0  sincb 

v N r * 'm  x Tx  m x Yx 
c 

(He.,e,  ) = x sin0„  - z cos0  cos<}> 

'■Neb  m x m xrx 

c 

(HeM^Q  ) = z cos0  sincj)  - y sin0 
v N 0 ' m x x 'm  x 

c 

3.3  Generation  of  Attitude  Data 


O.ice  the  axis  of  precession  is  determined,  then  the  generating  equation 
for  the  rocket  axis  X which  is  represented  by  e£  at  any  time  t,  is  given  by 
a slight  modification  of  equation  2 page  4 of  [4] . 

e^|  = ercosa  + evince  cos(wot+x)+c^sina(sinwot+x)  (3.3-1) 
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PP 


t,-1 

E 


masassBim, 


J 

A 


where 


3e 

Ga  = taq—  = -eQ  sin0  cosdi  -e , sin0  sind>  + e cos0 
0 30  ec  P rP  <J>C  P yp  rc  p 

V H5T  W - -\sl'%  * 


..  (H,eo) 

X = tan  1 Jr?-' — y"  -<»)  t at  a time,  t , of  the  maximum  of  axial 
v**e^J  0 0 o 

magnetometer 

rransforming  eJJ  to  the  North,  East,  and  vertical  system,  let 


a. ,=  e *eQ  a.,  = e *e 

11  r 0„  12  r d> 

c c. 


a21  " e0*e0  a22  = e0*e4> 
c c 


a ....  = e *e0  a--  = e »e 

al  <J>  0£  32  <j>  <j>c 


a,,  = e *e 
13  r r 


a23  = c0*er 


a„  ~ 

33  <J>  r 


then 


e£  = eQ  (a^cosa  + a2Jsina  cos  (u)Qt+x)  + a^sina  sin  fwot+x) 


+ 6,  (a^cosa  + a2,sina  cos  (a)Qt+x)  + a^sina  sin  (o3Qt+x)  (3.3-2) 

+ § (a^cosa  + a23sina  cos  (wot+x)) 
rc 


and  expressing  (3.3-2)  in  terms  of  direction  cosines 


where 


e"  = en  cosaQ  + e,  cosa,  + e cosa 
r 0 0 <f>  <f>  r r 

c c c 


cosa„  = e"*e 


(3.3-3) 


0 r 0. 


cosa.  = e”*6, 
<P  * <f> 


cosa  = e"*e 
r r r„ 
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can  also  be  expressed  in  terms  of  its  elevation  0a  and  azimuth,  <j>a 
angles  by 

e"  = §_  cos0  cos<j>  + e,  cos0  sin<f>  + e sin0  . 
r 0 a Ya  6 a ra  r a 
c Yc  c 

Then  by  equating  coefficients  of  (3.3-3)  with  (3.3-4),  the  attitude  at 
any  time,  t is  given  by 


sin0  = cosa 
a r 


(3.3-4) 


r1 

cosa. 


0 = sin-1 (cosa  ) 

a r 


tan<j>  = - 

Ya  cosa^ 

cosa, 

<p  = tan-1  ( - 

ra  vcosa„ 


3.4  Quality  Checks 

To  properly  check  for  the  best  determination  of  the  axis  of  precession, 
we  proceeded  in  the  following  manner.  First,  the  angles  between  the 
rocket  axis  and  the  Earth's  magnetic  field  vector  and  the  moon  were 
generated.  Then  the  times  of  the  minimums  and  the  maximums  of  the 
generated  data  were  found  and  the  comparison  of  these  times  with  the  ones 
from  the  actual  data  were  made  by  the  use  of  the  techniques  described 
in  Section  1.4.  The  axis  of  precession  which  generated  data  best  meeting 
the  experimenter's  criteria,  together  with  the  best  comparison  between 
the  actual  and  the  generated  outputs  was  selected. 


3.4.1  Generating  Magnetometer  and  Lunar  Data 

The  simulation  of  the  angles  of  the  rocket  axis  makes  with  the 
Earth's  magnetic  field  and  the  moon  was  performed  by  taking  scalar  products 
of  equation  (5.3-3)  with  (3.1-4)  and  (3.1-5),  i.e.. 


and 


H^(t) 


M”(t) 


Il*e"  = x cosa,.  + y cosa, 
r m 0 m q> 


+ z cosa 
m r 


M*e"  = m,cosan 
r 1 0 


m-cosa,  + m-cosa 
2 <p  3 r 


(5.4-1) 


(3.4-2) 
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3.4.2  Critical  Values  and  Second  Derivation  Test 

Once  the  data  is  generated,  the  times  of  the  maximum  have  to  be 
known.  First,  the  critical  values  H£(t)  and  M^J(t)  must  be  found.  This 
is  done  by  letting 

h = H*e 
r r 

Then  from  (3.3-1) 

h^(t)  = = h^cosa  + h^sina  cos  (wot+x)+h^sina(sintoot+x) 

and  the  critical  values  are  given  by 


he  * H'ee 


h.  » H*e, . 

4>  4> 


(3.4. 2-1) 


3h»(t) 


^ — = wosina[hQ-sin(a)ot+x)+h({)cos  (wQt+x]  = 0 


K -h 

i.e.,  when  tan(w  t+x)  = *—•  = -JP- 
0 n0  0 


From  [4]  by  the  second  derivative  test,  if  t is  the  time  when 

K 

ta"<“„Vx)  =5^-, 

9 


0 


* 


and 


+ h , 


sin  (w„t  +x)  = 


I L TV  1 - jmn  — ^ 

0 0 


9 <f>. 


then  t is  a time  of  a maximum  for  (3.4. 2-1) 


(3. 4. 2-2) 


(3. 4. 2-3) 


Substituting  (3.4. 2-2)  and  (3.4. 2-3)  in  (3.4. 2-1),  we  find  that 


h"(t  ) = h cosa  + sina\/h|  + h 
rv  o'  r >9 


0 ’ *> 


but 


implies 


h = cos$ 


H 


and 


4 


h2  -i  h?  + h2 
r 0 <p 


= 1 


4 


h9  * » * 


sinfr 


H 
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h"(t  ) = cos3„cosa  + sin3,,sina 
r o H H 

= cos(3h  - a). 


(3. 4. 2-4) 


Similarly,  if  is  the  time  when  h^(t)  is  a minimum  then 

h»(tj)  = cos(3h  + a).  (3. 4. 2-5) 

When  the  axis  of  precession  is  not  constant,  equations  (3.4. 2-4)  and  (3. 4. 2-5) 
do  not  hold  since  3jj  and  a then  become  functions  of  time. 

3.5  Modification  in  Analysis 

Using  the  preceeding  analysis  together  with  some  realistic 
lunar  sensor  angles  supplied  by  the  experimenter,  the  attitude  was 
determined  for  the  well-behaved  area.  Further  inspection  by  the  experimenter 
of  his  photometer  data  showed  that  the  initial  lunar  data  originally 
supplied  was  outside  the  acceptable  range  (-*4°) . From  the  structural 
pattern  of  the  photometer  data (Figure  24)  and  its  location  on  the  payload 
(Figure  l),it  was  possible  to  determine  an  approximate  angle  the  moon 
makes  with  the  lunar  sensor  at  certain  times.  These  determinations  all 
reflect  +4°  of  accuracy  and  +_.  5 seconds  in  time  due  to  the  variations 
in  the  telemetered  photometer  data.  The  existing  attitude  analysis 
was  modified  as  described  below  to  arrive  at  the  best  attitude  solution 
for  this  lunar  information. 

LetjT}  be  the  set  of  the  times  of  the  lunar  readouts  and  their 
associated  magnetometer  readings  in  one  second  interval  of  the  time 
supplied  from  the  photometer  in  formation. We  can  defincjT^as 

» • • • > • • • >tk[ 

i 

where  t ^ is  the  time  of  a lunar  readout  closest  to  the  time  supplied 
by  the  photometer  data. 
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and 


t1  - tpj.  - % second, 
tk  ~ tpk  + % second. 

Let  {l},  be  the  set  of  lunar  angles,  where  Lpk  is  the  approximate 
angular  value  derived  from  the  photometer  data  given  by 


where 


\l}  = 


Ll»...>Lpk>- 


4° 


Finally  let jTxL}  be  the  set  of  all  combinations  of  set  T and  L taking  two 
at  a time. 

{TxL[  - {(tijLj),  (tj  ,L2) , . . . , (tj , Lk) , . . .,(tk, Lk)  } • 

Then  the  attitude  was  calculated  for  each  of  the  points  of  the  setjTxl} 
by  equations  (2.1-7),  (2.1-8)  and  (2.1-9).  An  axis  of  precession  was 
calculated  for  each  of  these  orientations  by  (2.2-1)  and  (2.2-2)  and 
the  lunar  and  magnetic  field  information  simulated  by  (2.4-1). 

The  magnetic  pitch  angles  for  the  axial  magnetometerer,  generated 
by  using  each  of  the  calculated  axes  of  precession  discussed  above, 
were  compared  with  actual  magnetic  pitch  angle  (see  Section  1.4).  The 
orientation  of  the  axis  of  precession  that  produced  a minimum  P.MS  measure 
(sec  Section  1.3. 3.1),  that  satisfied  all  the  behavioral  patterns  of 
the  actual  magnetometer  data  and  all  the  experimenter’s  lunar  criteria, 
was  selected. 


-27- 


SECTION  4 


ATTITUDE  ANALYSIS  FOR  REMAINING  PORTION 
OF  THE  FLIGHT 


For  the  entire  time  the  photometers  functioned  properly  (i.e., 

51-376  seconds)  continuous  attitude  information  is  needed  even  when  the 
rocket's  axis  of  precession  is  not  constant.  This  requirement  means 
the  attitude  had  to  be  determined  for  certain  ascent  and  descent  portions 
of  the  flight  for  which  Section  3 does  not  apply.  The  approaches  to 
determine  the  attitude  differed  according  to  the  nature  of  the  attitude 
information  available. 

The  final  results  of  the  technique  described  below  in  addition  to 
those  results  previously  described  were  accepted  by  the  experimenter 
(sec  Figures  (25  - 30)). A particular  result  of  this  analysis, 

namely,  tne  angle  between  the  lunar  probe  and  the  moon  vector  compared 
favorably  (within  4°)  with  lunar  sensor  outputs  from  a similar  flight. 

4 . 1 Ascent  Portion 

For  the  ascent  portion  of  the  flight,  there  was  no  lunar  angular 
information  available.  To  determine  attitude  during  this  portion  of 
the  flight,  i.e.,  from  51  seconds  to  90.2  seconds  after  launch,  the 
orientation  of  the  rocket  axis  at  90.2  seconds  known  from  the  analysis 
for  the  well-behaved  portion  of  the  flight  was  used  as  a reference  position 
from  which  to  work  backwards.  The  apparent  precessional  motion  of  the 
axial  magnetometer  data  from  62.8  to  90.2  seconds  was  used  together  with 
(3.2-1)  and  (3.2-2)  to  determine  an  apparent  axis  of  precession  during 
this  period.  At  this  point,  we  could  simulate  lunar  and  magnetic  information 
with  respect  to  the  rocket  axis  by  using  (3.4-1)  and  (3.4-2)  respectively. 

In  a similar  manner,  i.e,,  using  the  generated  orientation  of  the  rocket 
axis  at  62.8  seconds  and  the  processional  motion  of  the  axial  magnetometer 
data  from  51  to  62.8  seconds,  the  attitude  was  determined  for  the  rocket 
axis  back  to  51  seconds.  The  results  in  each  of  the  intervals  discussed 
satisfied  the  data  requirements  of  the  axial  magnetometer  plus  all  the 
behavioral  characteristics  of  the  photometer  data. 
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4.2  Descent  Portion 

Unlike  the  ascent  portion,  two  types  of  lunar  angular  information 
were  available.  The  first  consisted  of  lunar  sensor  angular  outputs 
at  a few  discrete  points,  but  these  angles  reflect  a four-degree 
cone  of  uncertainty  (refer  to  Figure  31)).  The  attitude  was  then 
calculated  for  these  points  using  equations  (3.1-8)  and  (3.1-9).  The 
other  type  depended  upon  the  fact  that  the  axial  magnetometer  read 
approximately  zero  magnetic  field  at  about  360  seconds  after  launch. 
Then  equation  (3.1-7)  now  yields 

Therefore 


and 


Thus,  the  attitude  could  be  calculated  at  this  point  using  equation  (3.1-8) 
and  (3.1-9). 

Since  continuous  attitude  information  is  required  for  this  portion 
and  the  axis  of  precession  was  not  constant,  studies  showed  it  to  vary 
at  approximately  the  same  rate  as  the  midpoint  of  the  curve  of  the  axial 
magnetometer  output  (see  Figure  32)).  Then  using  this  information 
together  with  the  behavioral  pattern  of  photometer  data  and  the  discrete 
attitude  points,  the  techniques  discussed  in  Section  3 were  modified 
to  generate  the  attitude  information.  The  acceptance  tests  for  the 
attitude  data  here  were  the  same  as  for  the  well-behaved  area. 


(H.M)2  = Mullu  . 
(H*M), 


Mx  = sixi(cos_1(Mu)). 
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APPENDIX  A 


MAXIMUM  OUTPUTS  OF  THE  MAGNETOMETER  DATA 


A.  Specification  of  Curve  Fitted  Data  Areas 

Given  the  sinusoidal  nature  of  the  lateral  magnetometer  outputs, 
the  prepared  software  routines  have  the  capability  of  determining  either 
the  times  of  relative  maximum  outputs  or  the  times  of  relative  minimum 
outputs.  If  a least  squares  approximation  of  the  data  is  desired  in 
the  maximum  areas  of  the  output,  then  a test  is  performed  on  the  data 
satisfying  the  inequality 

B (T . ) < V. . 
i'  1 

Likewise,  if  a fit  of  the  minimum  areas  is  desired,  then  only  those 
values  satisfying  the  inequality 

Bfr.)  > V. 

are  stored  for  computation  in  the  least  squares  polynomial  fits  of  the 
specified  areas  of  the  data. 

B.  Acceptance  Level  for  the  Least  Squares  Data 

In  order  to  initiate  the  procedure  for  storing  data  in  the 

maximum  areas  of  the  least  squares  approximation,  a value  V,.  must  be 

found  such  that  V.  > Bft.l  and  V.  .<  B(t.  Once  initiated,  data 

i v i'  i-l  — i-l 

is  stored  until  V.  < B(t.)  and  V.  . > B(t.  .)  where  j > i.  This  last 
1 J 1-1  - 1-1 

condition  is  a necessary  one  for  ending  the  sample  region  but  not  a 
sufficient  condition.  The  sufficiency  test  is  made  by  examining 
V.  . to  ensure  that  V.  is  not  a noise  point.  Tf  V.  . < B(t.  .)  then 

j+1  j 1 j + 1 j+1 

a least  squares  parabolic  fit  is  performed  on  the  sampled  data.  However, 
if  Vj+1  B(tj+1),  the  point  V.  is  discarded  and  the  summing  procedure 
of  the  least  squares  continues.  Similar  tests  are  made  in  the  minimum 
regions  of  fitted  data  except  that  the  inequality  signs  arc  reversed. 

C.  Least  Squares  Approximation 
1.  1st  derivative  test. 

In  each  specified  region,  the  data  is  fit  to  the  quadratic 


FT 
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The  value  of  the  bias  B(t^)  at  which  the  sampling  begins  and  that  value 
B(tj),  j > i at  which  the  sampling  ends  are  obviously  just  the  points  of 
inflection  for  the  sinusoidal  curve.  The  time  of  relative  maximum 
(minimum)  is  obtained  by  setting  = 0 and  solving  for  time  t.  This, 
in  turn,  provides  us  with  a predicted  output  value 

Y = At£  + BT„  + C 
M M 


where 


T = 


B ' 


2K 


2.  Error  tolerance 


Should  the  measured  and  predicted  outputs  deviate  by  more  than  some 
specified  limit,  usually  taken  to  be  4%,  the  acceptance  level  for 
least  squares  data  may  be  modified  so  that  the  sampled  regions  more 
closely  approximate  the  parabolic  fit.  This  modification  amounts  to 
accepting  the  area  containing  maximum  data  when  > B(t^)  + k where  k 
is  a constant  determined  from  examination  of  the  oscillographs.  For 
minimum  output  regions,  the  inequality  takes  the  form  V.  < B(t^)  - k. 


3.  Output  quantities 

For  each  value  of  T^  generated  by  the  1st  derivative  test,  an 
associated  Y(T^)  value  and  AT^  are  output.  Special  note  is  made  of  any 
value  of  AT^  outside  a specified  interval  about  the  assumed  rotation 
rate  and  further  investigations  of  these  areas  are  undertaken  as 
necessary.  A mean  AT^  value  is  computed  for  the  well-behaved  portion 
of  the  flight  and  output  at  the  end  of  the  data  pass.  This  AT^  value 
for  each  lateral  magnetometer  normally  consists  of  spin  and  process jonal 
motion  and,  therefore,  is  referred  to  as  the  rotation  rate. 


Although  no  specific  mention  as  yet  has  been  made  of  the  axial 
magnetometer  output,  the  routines  apply  equally  as  well  in  the  case  of 
this  magnetometer.  Computation  of  maximum  and  minimum  values  provide 
, us  with  a good  approximation  to  the  cone  angle  of  the  vehicle  axis. 

t 

■ Should  the  output  exhibit  only  the  slightest  variance  in  successive 

| maximum  and  minimum  values  respectively,  then  the  vehicle  axis  motion 

• is  due  almost  entirely  to  precessional  motion  and  the  cone  angle  will 

remain  nearly  constant. 
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If  the  vehicle  axis  study  option  is  desired  in  addition  to  the  lateral 
magnetometers,  supplemental  outputs  include  the  cone  angles  and  precession 
period  of  the  vehicle  axis. 
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INSTRUMENT  ORIENTATIONS 


. LUNAR  SENSOR 


^ 9 (Azimuth) 

V 


(View  from  nose) 


AZIMUTH 

INSTRUMENT (degrees) 

X Magnetometer  0 

Y Magnetometer  225° 

Z Magnetometer  315° 

Lunar  Sensor  0 

2750  A Photometer  0 

2600  A Photometer  0 

Siemens  Brushless  Motor  90° 


ELEVATION 

(degrees) 

0 

90° 

90° 

90° 

85° 

80° 

30° 


Figure  1 
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COORDINATE  SYSTEM 
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Figure  2 
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Figure  3 
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350 


ANGLE  OF  ATTACK 


0.0  150.0  200.0  250.0  300.0  350.0  400.0 

SEC0NOS  AFTER  LAUNCH 


RAM-5C 


SPECIFICATIONS 


Figure  18 


Input  voltage:  24  to  32  volts  dc 
Input  current:  11  milliarnpcres  (typical) 

Range  of  field:  0 to  ±600  milliocrsteds 
Sensitivity:  .004  volts  dc  per  millioersted 
Stability  of  sensitivity:  ±3% 

DC  output  for  zero  field:  2.40  ±.02  volts  dc* 

Stability  of  output  for  zero  field:  ±.025  volts  dc 

Operating  temperature  range:  165°F  to  0°F  (±185°F  to  — 50“F  vich  reduced  accuracy) 

Linearity:  ±3%  of  full  scale 
Frequency  response:  195  Hz  (typical) 

Output  impedance:  less  than  20  kfl 
Output  load  (for  factory  calibration):  100  kfl 
Length  of  senor  cable:  5 feet 
Weight:  5 ounces 

•Standard  units  are  biased  as  specified.  An  instrument,  designated  RAM-5C  NB,  is  available  with 
factory  modification  to  provide  0-volt  output  for  zero  field. 


SCHONSTEDT  INSTRUMENT  COMPANY 

1775  WIEHLE  AVENUE 
RESTON,  VIRGINIA  22070 
Area  Code  703  471-1050 
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Figure  19 


.CALIBRATION  nATA 
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Output  Signal 
in  Volts  DC 

v.flo 

V-  bA 

V.V 3 

V.23 

-JLQjL 

3.fl/ 

3.  VO 


-^3 ,10 

3.00 


2, VO 


-.7? 


Heli  Flux' 

MAGNETIC  ASPECT  SENSOR 
TYPE  RAM- SC 


(Bias  Level) 


Direction  of  Magnetic  Field  for 
Voltage  Signals  Above  Bias  Level 


NOTE: 

Calibration  made  with  a 100k 
OHM  resistor  from  signal 
output  to  negative  terminal 
of  battery  source,  and  a 100k 
OHM  resistor  from  bias  output 
to  negative  terminal  of  battery 
source. 


SCHONSTEDT  INSTRUMENT  COMPANY 
RESTON,  VIRGINIA 
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LSQMAG  FLOWCHART 


START 


RMS  . 

= 0 

mm 

k 

= 0 

Initialize  cx^  i=l,2,3 
Iteration  limit  L,  increment  value  m 


Compute  the  least  squares  terms  in  the  system  established  by 

31 


where 


In  - 2[H2-A1(XBg.aJ)!  - A2(YB]g*a2)2  - A3(A 


Compute  RMS  =^f — - Xi 


where  X^  is  the  difference  between  measured  and  theoretical  field 

/ 

1 — ■ - r?\ 


Is  RMS  . < RMS 


Update  the  a., 
and  + Va7  values 


Is  k > 2L 
YES 

_J 

CALL  EXIT 


Figure  23 


A09. 209-1  LUNAR  ANGULAR  READOUTS 


L 


* 

■ 


1 

J 
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sec. 

+10° 

11 
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min. 
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+10° 

11 

hrs . 

10 

min. 

27.05 

sec. 

+ 5° 

11 

hrs. 

10 

min. 

27.35 
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0° 

Note:  these  outputs  are  all  within  + 4°  of  the  true  angular  value 


Figure  29 


A09.209- I 

AXIAL  MAGNETOMETER  CORRECTED  OUTPUT 
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